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SCALING LIMIT OF FLUCTUATIONS IN STOCHASTIC 
HOMOGENIZATION 
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Abstract. We investigate the global fluctuations of solutions to elliptic equa¬ 
tions with random coefficients in the discrete setting. In dimension d ^ 3 and 
for i.i.d. coefficients, we show that after a suitable scaling, these fluctuations 
converge to a Gaussian field that locally resembles a (generalized) Gaussian 
free field. The paper begins with a heuristic derivation of the result, which can 
be read independently and was obtained jointly with Scott Armstrong. 
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1. Heuristics 

The goal of this paper is to give a precise description of the fluctuations of 
solutions of elliptic equations with random coefficients, in the large scale limit. 
Before stating our precise assumptions and results, we present powerful heuristics 
that enable to guess the results and give a better comprehension of the phenomena 1 . 
These heuristics were obtained in collaboration with Scott Armstrong, whom we 
warmly thank for letting us include this material here. 

1.1. The (generalized) Gaussian free field. We start by introducing white noise 
and (generalized) Gaussian free fields. These will be the fundamental objects used 
in the heuristic derivation of the large-scale behavior of the first-order correction to 
stochastic homogenization below. 

The random distribution w is a (one-dimensional) white noise with variance a 2 if 
for every (j) 6 (^(K), «;(</>) is a centered Gaussian random variable with variance 
a 2 f (f> 2 . (We can in fact define w(4>) for any <f) 6 L 2 (R d ) by density.) Informally, 

E[u>(x) w(y )] = cr 2 5{x - y), 

where S is a Dirac mass at the origin. More generally, the random, d-dimensional 
distribution W = (Wi W,i) is a white noise with covariance matrix Q if for every 

</>=(^,...,</> d )6C~(R d ), 

W(<t>) -WM + - + WM 

is a centered Gaussian random variable with variance f (f> ■ Q(/). Informally, 

E [Wi(x)Wj{y)] = QijS(x-y). 

In dimension d = 1, one way to define a Brownian motion B is to ask it to satisfy 
(1.1) B' = w, 

where w is a one-dimensional white noise, and B' denotes the derivative of B. The 
Gaussian free field is a high-dimensional version of Brownian motion. By analogy 
with (1.1), we may want to ask a Gaussian free field <1? to satisfy V*! 1 = W, where W 
is a d-dimensional white noise. However, this does not make sense because W is 

1 The recording of a talk presenting this is also available at http://goo.gl/5bgfpR. 
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not a gradient field; so we will instead define V4 1 as the L 2 projection of W onto 
the space of gradient fields. In view of the Helmholtz-Hodge decomposition of any 
vector field into a potential part and a solenoidal part: 

(1.2) L 2 = (Vrt) © {g : V • g = 0}, 
this leads us to the equation 

(1.3) -a $ = v-w: 

A minor variant of the construction above is to consider the Helmholtz-Hodge 
projection with respect to a uniform background metric given by a symmetric 
matrix Oh- In this case, equation (1.3) becomes 

(1.4) - V-a h V4> = V-W. 

Since the equation is linear, one can give a mathematically precise sense of V$ using 
classical arguments of the theory of distributions. (In dimension one, this also gives 
a sensible way to define a “Brownian motion” on the circle, i.e. a Brownian bridge, 
as opposed to (1.1) which does not satisfy the compatibility condition f w - 0.) We 
take (1.4) as the definition of the Gaussian free field associated with ah and Q. 

On the full space with d ^ 3, one can also define $ itself (and not only V4?), 
e.g. as the limit as /i tends to 0 of <J> M such that 

(At- V-a h V)3v = V-W. 

In this case, one can express the two-point correlation function of in terms of the 
Green function Q h of -V • cthV and the covariance matrix Q of W: 

(1.5) E[4>(0 )$(x)]=J vGh(y)-QvGh(y-x)dy. 

If Q happens to be a multiple of Oh, then an integration by parts enables to replace 
the integral above by a constant times Gh(x), and we recover the more common 
definition of the Gaussian free field as a Gaussian field whose covariance kernel is 
a Green function. However, for generic ah and Q, the correlation in (1.5) cannot 
be expressed as a Green function. In other words, our definition of (generalized) 
Gaussian free field is wider than the standard one. 

It is important for the remainder of the discussion to be familiar with the scaling 
and regularity properties of white noise and Gaussian free fields. As for white noise, 
W(r ■) has the same law as r~^W. In particular, thinking of r 0, we see that 
zooming in on W at scale r produces a blow-up of (and conversely if we think of 
r -s- +oo). This is an indication of the fact that W has (negative) Holder regularity 
a for every a < -f, and no more. In view of (1.4), the Gaussian free field <f> is 
such that 4>(r •) has the same law as r~^ +1 4>, and has Holder regularity a for 
every a < — | + 1. In dimension d = 1, we recover the fact that Brownian motion 
trajectories have Holder regularity a for every a < ^. In higher dimensions, the 
Gaussian free field fails to have regularity 0; it only makes sense as a distribution, 
but not as a function. 

1.2. Homogenization and random fluctuations. We now turn to the homoge¬ 
nization of the operator -V • aV, where a : R d R dx<i is a random field of symmetric 
matrices. We assume that the law of a is stationary and posesses very strong mixing 
properties (e.g. finite range of dependence), and that Id ^ a ^ Cld for some constant 
C < oo. In this case, it is well-known that the large scale properties of the operator 
-V • aV resemble those of the homogeneous operator -V • V, for some constant, 
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deterministic matrix a h. Our goal is to describe the next-order correction. For 
p € R d and := |) d £ R d , we introduce 


( 1 . 6 ) 


v{n r 


,p) := inf 
K (Dr) 


L 


-(p+ Vu) 


a(p + Vu) 


(where / n = |□ r .| 1 f D ). This quantity is subadditive: if is partitioned into 
subcubes {y + □ ,,.),,, then u(Q r ,p) is smaller than the average over y of v(yy + IH s ,p). 
Indeed we can glue the minimizers of v(y + IH s ,p) and create a minimizer candidate 
for u(D r ,p). Roughly speaking, it was shown in [6] that homogenization follows 
from the fact that 


(1.7) 


u(Q r ,p) -► ~P' a bP 

r-*- oo 2 


(which itself is a consequence of the subadditve ergodic theorem). It is natural to 
expect the next-order correction to homogenization to follow from the understanding 
of the next-order correction to (1.7). However, the next-order correction to (1.7) is 
driven by a boundary layer, which is of order r -1 (see [2]), and is not relevant to the 
understanding of the interior behavior of solutions. We thus assume that v has been 
suitably modified into v in order to get rid of the boundary layer. After performing 
this modification, we expect u( • ,p) to be “almost additive” [2], and therefore that 

\U r \ 112 [v(U r ,p) - -p-a h pj 


converges to a Gaussian random variable as r tends to infinity (as a consequence of 
the strong mixing assumption on the coefficients). Closely related statements were 
proved in [24, 4, 27, 25, 11]. We want to encode this information in a way that is 
consistent with respect to changing the vector p, the scale r and translations of the 
cube. For this purpose, we let IF be a matrix-valued white noise field such that as 
r becomes large, 

(1.8) v(x + U r ,p) =: -p- (a h + W r (x))p, 
where W r is the spatial average of IF on scale r: 

(1.9) W r {x) := / IF. 

J x+n r 

(We may also think of W r as the convolution of IF with a rescaled bump function: 
Wr := IF * x (r) , with X 6 C c °°(R d ,R + ) such that / % = 1 and x (r) == r~ d X (-/r).) This 
encodes in particular the fact that v(x + □ r ,p ) and v(yy + O r ,q) are asymptotically 
independent if x + and y + n r are disjoint. Recall that each coordinate of W r (x ) 
is of order |n r | — = r~ d ! 2 « 1. We interpret (1.8) as indicating that if a function 
locally minimizes the energy over x + D r and has average gradient p, then its energy 
over x + is approximately \p- (ah + W r (x))p. 

The corrector for -V • aV in the direction p is usually defined as the sublinear 
function solving 

-V • a(p + V<£) = 0 

in the whole space. We think of <f> as the minimizer in the definition of u(D/?,p), 
for R extremely large (in fact, infinite), and focus our attention on understanding 
the spatial average (f> r of </> on scale r, 1 « r « R. The discussion above suggests 
that (f) r minimizes the coarsened energy function 

\{P+ Vv) • (a h + W r )(p + Vu), 

JDr 2 

whose Euler-Lagrange equation is 


-V • (ah + W r )(p + V<(v) = 0. 
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Rearranging, we obtain 

-V • (a h + W r )S/4> r = V • (W r p). 

Since W r « 1, we have V</v « 1- Therefore, the term W r on the left-hand side can 
be neglected, and we obtain the Gaussian-free-field equation 

(1.10) - V-a h V</v = V-(W r p). 

We summarize this heuristic computation as follows. 

• We let W be the matrix-valued white noise field whose covariance is related 
to the fluctuations of the energy via (1.8); 

• we let $ be the random distribution defined by 

(1.11) - V-ahV$ = \/-iWp) ; 

• then the large-scale spatial averages of the corrector </> have about the same 
law as those of <f>. In other words (and by the scale invariance of $), the 
random distribution r 5_1 <)(r •) converges in law to $ in a suitably weak 
topology, as r tends to infinity. 

A similar analysis can be performed for, say, solutions of equations of the form 

(1.12) - V • aVu = / in R d , d ^ 3, 

where / € C™ (R d ) varies on scale £ _1 » 1. (The function / should be of order e 2 in 
order for u to be of order 1.) We consider the spatial average u r of u over scale r, 
1 « r « s~ 1 . By the same reasoning as above, we expect u r to satisfy the coarsened 
equation 

-V • (a h + W r )S7u r = f. 

We write u r = uu + u r , where u h solves 

-V-a h Vuh = /, 

so that 

-V • (tth + W r )\/u r = V • (W r V«h)- 

As before, we expect the term W r on the left-hand side to be negligible, so we obtain 

(1.13) - V-OhVu r = V-(W r Vw h ). 

We stress that if instead we use the formal two-scale expansion u - + T,i (p^^iUh 

and the large-scale description of the corrector in (1.10), we are led to a different 
and incorrect result. 

We can again summarize our conclusions as follows. 

(1) We define the matrix-valued white noise field as above, according to (1.8); 

(2) we let be the random distribution defined by the equation 

(1.14) - V • a h V^ = V • (ITVuh) ; 

(3) Then the large-scale averages of u - have about the same law as those 
of 

The random distribution <2/ is not a (generalized) Gaussian free field per se because 
the term Vrth appearing in its defining equation varies over large scales. However, 
it will have similar small-scale features. If we normalize / so that u h is of order 1, 
then Vuh is of order e, over a length scale of order e -1 . Hence, we can think of 
as locally like a (generalized) Gaussian free field (of order 1), and being close 
to 0 outside of a domain of diameter of order e -1 . Incidentally, if we make the 
(unjustified) ansatz that u - uy, is approximately a regularization of % on the unit 
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scale (that is, if we pretend that the conclusion (3) above actually holds on the unit 
scale), then we recover the (correct) error estimate 


f |w-wh| 2 j 


1/2 


elog 1/2 (e: _1 ) ifd = 2, 
£ if d ^ 3. 


Mathematically speaking, a version of the statement that r2 -1 0(r •) converges 
to a (generalized) Gaussian free field was proved in [23, 22], The aim of the present 
work is to prove a version of the conclusion (3) above. 


Before turning to this, we want to emphasize why we believe these results to be 
of practical interest. Homogenization itself is interesting since it enables to describe 
approximations of solutions of equations with rapidly oscillating coefficients by 
solutions of simple equations described by a few effective parameters. What the 
above arguments show is that the same is true of the next-order correction: in order 
to describe it, it suffices to know the few parameters describing the covariance of 
the white noise W. The white noise W takes values in symmetric matrices, so its 
covariance matrix is fully described by N(N + 1)/2 parameters, where N = d(d+ 1)/2 
(6 parameters in dimension 2, 21 parameters in dimension 3). Naturally, fewer 
parameters are necessary for problems with additional symmetries. For instance, the 
noise is described by only one parameter in any dimension if we assume rotational 
invariance. 


2. Introduction 

2.1. Main result. We focus on dimension 3 and higher and on a discrete setting. 
Our main assumption is that the random coefficients are i.i.d. and bounded away 
from 0 and infinity. Our goal is to justify the points (1-3) listed above in this 
context. 

In order to state our assumptions and results more precisely, we introduce some 
notations. We work on the graph (Z d ,B) with d ^ 3, where B is the set of nearest- 
neighbor edges. Let e \,..., be the canonical basis in Z d . For every edge e e B, 
there exists a unique pair (e, i) e Z d x {1,..., d} such that e links e to e + e^. We 
will write e = e + and e = (e, e). 

We give ourselves a family of i.i.d. random variables indexed by the edges of the 
graph. For convenience, we will assume that these random variables can be built 
from a family of i.i.d. standard Gaussians. More precisely, we let £ = (Ce)edB 6 H := R B 
be i.i.d. standard Gaussian random variables, and rj be a smooth function that is 
bounded away from zero and infinity with bounded first and second derivatives. 
The conductance associated with the edge e is then defined to be ?/(Ce)- The space 

is equipped with the product Borel cx-algebra, and we denote the law of £ by P, 
with associated expectation E. 

Let a : Z d x 12 R dx “ be the matrix-valued function such that a(x,() = 
diag(ai(a:, £),..., dd(x, £)) with di(x,() - v(((x,x+et ))■ I* 1 what follows, we will 
most of the time keep the dependence on £ implicit in the notation. 

For any / : Z d -» R we define the discrete gradient V/ := (Vi/,..., Vd/) by 
Vi/( x) = f(x + 6i) - f(x). For any g : Z d -»■ we define the discrete divergence 
:= Eli by V*gi(x) = g t (x - e*) - gi(x). We define V e , V* similarly for 
functions defined on eZ d , i.e., x) = e~ 1 (h(x + se i ) - h(x)) and V* ih(x) = 

£^ 1 (h(x - set) - h(x)) for h : eZ d -*■ R. 

For any £ > 0, we consider the following elliptic equation with a slowly varying 
source term: 

(2.1) 


\/*a{x)SJu{x) = f{ex) (xeZ d ), 
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where / e C c (K d ) is compactly supported and continuous. The unique solution to 

(2.1) that decays to zero at infinity is given by 

u(x) = Y, G(x,y)f(ey), 

yeZj d 

with G(x,y ) the Green function of V*a(x)V (recall that the dependence on f is 
kept implicit, so u and G are random). 

We define u E (x) = £ 2 u(|), which solves 

V*a(-)V £ u e (x) = f(x) (xeeZ d ), 

£ 

and its limit u h, the solution of the homogenized equation in continuous space: 

-V • OhVuh(x) = f(x) (xeK d ), 

where the homogenized matrix a h is deterministic and constant in space. 

We are interested in the random fluctuations of u e after a spatial average. In 
other words, we think of u E as a (random) distribution U E J , which acts on a test 
function g e C c (R d ) as 

(2.2) U { J\g):=£ d £ u E (x)g(x) = e d+2 £ G(x,y)g(ex)f(ey). 

xeeZj d x,yeZj d 

Following [21], for any a < 0, we denote by = Cj“ c (R d ) the (separable) local 
Holder space of regularity a. Here is our main result. 

Theorem 2.1. Recall that d ^ 3. We define the random distribution ' E ^ by 

%U\g)=e-i{Ug\g)-nW ) {g)}). 

For every a < 1 - the distribution ^/ E ^ converges in law to as e -»■ 0 for the 

topology o/Cj“ c , where is the Gaussian random field such that for geC c (R d ), 

<% u \g) is a centered Gaussian with variance 

(2.3) a 2 g = JG f (x, z)g(x)g(z)dxdz 
with JGf given explicitly by (3.7). 

Remark 2.2. We can write the solution to (1.14) formally as 
&(x) = f Gh(x-y)V ■ (W(y)\7u h (y))dy 

jR d 

= f VGu(x - y) -W(y)\7Gh(y ~ z)f(z)dzdy, 

jR 2d 

where G h is the Green function associated with -V • ahV. By a straightforward 
calculation, this random distribution tested with g has a Gaussian distribution with 
mean zero and variance given by 

E{| f (x)g(x)dx\ 2 } = f JFr(x,z)g(x)g(z)dxdz, 

JR d JR 2d 

provided the covariance of the white noise W is given by {Kj^j } that appears in 
the definition of Jfj(x,z). In other words, the limiting distribution obtained 

in Theorem 2.1 can be represented as the solution of (1.14). 

Remark 2.3. Theorem 2.1 implies the joint convergence in law of (fiV/Y' ■*,... 
to a Gaussian vector field whose covariance structure is obtained by polarization of 
f i-*. jFf. Indeed, this follows from the fact that Theorem 2.1 gives the scaling limit 
of any linear combination of 9/^' ■*...., by linearity of / i-> a i/ E ^ ■ 
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Remark 2.4. A similar result is proved in [15] when d = 1, using a different method. 
In this case, ^ in (1.14) has Holder regularity a for every a < so it makes sense 
as a function. We expect the result to hold for d = 2 as well, but our method would 
have to be modified to handle the fact that only the gradient of u is really defined 
by (2.1) in this case. 

Theorem 2.1 is a consequence of the following two propositions. 

Proposition 2.5. For every g 6 C c (R d ), ^^\g) converges in law to ^/^\g) as 
e ->■ 0 . 

Remark 2.6. In fact, when a 2 + 0, our proof gives a rate of convergence of ty^\g) 
to ?/ (/) (s), see Remark 4.8 below. 

Proposition 2.7. For every a < 1 - is tight in C j“ c . 

2.2. Context. Stochastic homogenization of divergence form operators started 
from the work of Kozlov [17] and Papanicolaou-Varadhan [26], where a qualitative 
convergence of heterogeneous random operators to homogeneous deterministic ones 
is proved. Quantitative aspects were explored as early as in [28]. However, optimal 
bounds on the size of errors were obtained only recently in a series of papers 
[12, 13, 20, 8, 9, 14, 7[. Regularity estimates that are optimal in terms of stochastic 
integrability have been worked out in [1, 3, 9[. 

Our focus in this paper is to go beyond estimating the size of the errors, and 
understand the probability law of the rescaled random fluctuations. In this direction, 
central limit theorems for approximations of homogenized coefficients are obtained 
in [24, 4, 27, 25, 11]. The scaling limit of the corrector is investigated in [23, 22]. In 
the continuous setting, [16] indicates that when d ^ 3, the corrector should capture 
the first order fluctuation of the heterogeneous solution in a pointwise sense, but 
it is not clear whether it captures the fluctuations of the solution after a spatial 
average. A surprising feature of our result is that the limiting fluctuations are not 
those induced by the corrector alone. 

Our approach is based on that of [23, 22]. The fact that '(g) divided by its 
standard deviation converges to a standard Gaussian is derived using a second order 
Poincare inequality developed by Chatterjee [5], in the spirit of Stein’s method. (We 
will in fact use a slightly more convenient form of this result derived in [22].) The 
main difficulty lies in the proof of the convergence of the variance of ?/) '(g). A 
Helffer-Sjostrand formula enables to rewrite this variance in terms of gradients of the 
Green function. A quantitative two-scale expansion for the gradient of the Green 
function was worked out in [23]. Here, we follow the idea of [10] of introducing 
a stationary skew-symmetric tensor, which is denoted by {<Jijk}i j k=i an d relates 
to the flux (see Lemma 4.5). In the language of differential forms, the flux in the 
*-th direction is a co-closed 1-form, and we represent it as the co-differential of the 
2-form crj. This object enables us to represent the error in the two-scale expansion in 
divergence form, and thus significantly improve the two-scale expansion and simplify 
the subsequent analysis. 

2.3. Organization of the paper. The rest of the paper is organized as follows. We 
introduce basic notation and recall key estimates on correctors and Green functions 
in Section 3. Then we present some key ingredients in proving Theorem 2.1 in 
Section 4, including the Helffer-Sjostrand covariance representation, a quantitative 
two-scale expansion of the Green function and a second order Poincare inequality. 
The proofs of Propositions 2.5 and 2.7 are contained in Sections 5, 6 and 7. Technical 
lemmas are left in the appendix. 
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3. Setup 

3.1. Asymptotic variance. For x 6 Z d , we define the shift operator t x on S! by 
(t x C) e = Cx+ei where x + e := (x + e, x + e) is the edge obtained by shifting e by x. 
Since {/ e } e eB are i.i.d., {r x } x€ id is a group of measure-preserving transformations. 
With any measurable function / : VL -* R, we can associate a stationary random field 
f(x,() defined by 

(3-1) />,C) = (T x /)(C) = /(r x C). 

The generators of T x , denoted by {Th}f =1 , are defined by Dif := T ei f-f. The adjoint 
D* is defined by D* f := T_ e .f-f. We denote the gradient on by D = (D 1; ..., Dj) 
and the divergence D*g := T,f=i D*gi for g : fl -* M. d . The inner product in L 2 (0) 
and norm in L P (H) are denoted by (•,•) and || ■ || p respectively. 

Most of the time, we keep the dependence on / implicit and write f(x) = = 

f( T xO- For any e € B, the discrete derivative on e is defined by V/(e) := /(e) - /(e). 
If e is in the z-th direction, we define /(e) := /, for any / € K d , as the projection of 
/ onto e. 

For the random coefficients appearing in (2.1), we can define 

a(C) = diag(«i(C),...,a d (C)) := diag(r?(C ei ),... ,r?(Ce d )) 

so that a(x,() = cl(t x () = diag(? 7 (/ a:+ei ),... ,r](( x+ed )). Note that we also used 
ei,...,ed to denote the corresponding edges (0, ei),..., (0, e^)- Recall that we 
assume that (W 1 < g < C and \g'\, \q"\ < C for some C < oo. For simplicity we will 
henceforth write a e = g(Ce)- 

Under the above assumptions, it is well-known that there exists a constant matrix 
ah such that the operator V*aV homogenizes over large scale to the continuous 
operator -V • OhV, the Green function of which we denote as 

One of the main ingredients in the analysis of stochastic homogenization is the 
so-called corrector. For any fixed / 6 M. d and A > 0, the regularized corrector <j)\^ is 
defined through the following equation on probability space: 

(3.2) + D*a(D(f)\£ + /) = 0. 

It is proved in [12] that as A -»■ 0, <f>\^ -* <j>^ in i 2 (U), i.e., a stationary corrector <j>£ 
exists such that 

(3.3) D*a(D(/) ( + 0 = 0. 

For i = 1,..., d, we will write (f>i = <f> ei and = </> x, ei - The homogenized matrix ah 

is given by 

(3.4) / T a h / = E {(/ + H/» 4 ) r a(/ + H^)}. 

In the context of i.i.d. randomness, we have eth = aid for some constant a, where Id 
is the identity matrix. 

For a random variable F e L 2 (U), we say that U = d e F 6 L 2 (H) is the weak 
derivative of F with respect to if the following holds: for any finite subset A € B 
and any smooth, compactly supported function G : -»■ K, we have 

r)G 

(3.5) E{U G(/)} = E{FCeG(C)} - E{F—(/)}, 

d(e 

where G(/) depends only on {Ce'je'eA- We also call d e the vertical derivative , and 
by (3.5), its adjoint (under the Gaussian measure) is 


dt = ~d e + Ce- 
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We define df := (d e f) es B and for F = (F e ) eiM , d* F := £ esB 9*F e . The vertical 
Laplacian on the probability space is then defined by 

JSf := d*d. 

For i,j,k,l = 1,..., d, define 

Kijki(e) := (d e a e (ei+V<^i)(e)(ej+V</>j)(e), (l+.if)“ 1 <9 e a e (efc+V0fc)( e )( e / + V^i)(e)), 
and 

d 

(3-6) Kijki := X Kij k i(e n ). 

n= 1 

The kernel Jtf(x,z) appearing in Theorem 2.1 is given by 

(3.7) 

d z' 

Jff(x,z) = X Kijki d Xi g h (v - x)d Xo Gu(v - y)d Xk g b (v - z)d xi g h (v - w) 

i,j,k,i =l 

f(y)f (w)dydwdv. 

3.2. A comparison with the two-scale expansion. Let us see that the global 
fluctuations of w £ are not those suggested by its two-scale expansion. Recall that u e 
and u b satisfy 

V* £ a(-)V £ u £ (x) = f(x ) (x e eZ d ), 

£ 

and 

-V • a b Vuh(x) = f(x) (xeR^). 

A formal two-scale expansion gives 

X 

(3.8) u e (x) = u b (x) + eVuh(x) •</>( — ) + o(e), 

E 

where </> = (0i,..., (f>d)- In the continuous setting, (3.8) is proved rigorously with 
o(e)/e -» 0 in L l { fl ) for fixed x [16, Theorem 2.3], i.e., the first order correction is 
indeed given by the corrector in a pointwise sense. Since <f> is centered, we have a 
similar expansion for the random fluctuation, i.e., 

X 

u e (x) = E{u e (x)} + eVwh(^) • 0( —) + o(e). 

£ 

Concerning global fluctuations, we need to compare the random field u e (x) - 
E{it e (:r)} with s\Ju b (x) ■ 0(f)- For a test function g 6 C c (R d ), Theorem 2.1 shows 

(3.9) £5 XI ( u e(x) -E{u e (a;)})g(x) => N(0,a 2 g ), 

xzelA 

with <jg = / R2 d o#/ (x, z)g(x)g(z)dxdz. By [22, Theorem 1.1], we have 
£5 X eVuh(x) -<t>(-)g(x) => N(0,a%), 

with of = / R2 d J£/(x, z)g(x)g(z)dxdz and 

(3.10) 

d r* 

Kijki I &xi gh (v ~ x^dxj g b (x — y') d X k gh (v ~ z')d xi g b (z — w') 
i,j,k,l=1 JR 

f(y)f (w)dydwdv. 

If of and of were equal for every admissible / and g , then 
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and 

d r* 

'y ^ R-ijkl I ^Xi Gh (y ~ X^9 X j Gh (*£ — y)^:Efc Gh — Gh (% ~ dv 

i,j,k,l= 1 jRi 

would have to be equal almost everywhere (as functions of x^y,z and w ). However, 
the first quantity diverges when y gets close to w since d ^ 3, while this is not so for 
the second quantity. This shows that the fluctuations of u e are not those suggested 
by the two-scale expansion. 

The heuristics in Section 1 provide a clear picture of the above phenomenon 
(other than the explanation that o(e) may contribute on the level of £5 when d ^ 3). 
If we write the solution to (1.11) as <b(;r) = f Rd \7Gh(x - y) ■ W(y)pdy, the rescaled 
limit of the corrector eVMh(^) ■ <(>(f) is 

(3.11) £~z +1 Vu h (x)f VGh(x-y) • W(y)Vuu(x)dy, 

£ J M d 

and we already know from Remark 2.2 that the rescaled limit of u e (x) -E{it £ (x)} is 

(3.12) £~%(u e (x) -E{u £ (x)}) f VGh(x - y) ■ W(y)Vu h (y)dy. 

J R d 

Comparing the r.h.s. of (3.11) and (3.12), it is clear that they are two different 
Gaussian random fields. We further observe that they are linked through a Taylor 
expansion of Vwh(y) around x. By writing 

Vu h (y) = Vri h (z) + V 2 u h (x)(y-x) + ..., 

with V 2 Uh(."c) the Hessian of Uh, we have 

f VGb(x-y) ■ W(y)Vu h (y)dy 

jR d 

(3.13) = f VGh(x- y) ■ W(y)\/Ub(x)dy 

jR d 

+ [' VGh(x - y) -W(y)V 2 Uh(x)(y - x)dy + .... 

jR d 

The term f Rd VGh(x - y) ■ W{y)'\J 2 u^{x){y - x)dy should correspond to the second 
order corrector obtained by the two-scale expansion, and we also expect those higher 
order terms appearing in (3.13) to correspond to the rescaled limit of the higher 
order correctors (provided that they are stationary). It does not seem possible to 
simply add a finite number of terms in the two-scale expansion to recover the correct 
limiting field. 

3.3. Properties of correctors and Green functions. We summarize here sev¬ 
eral results obtained in [12, 18] which will be used frequently throughout the paper. 
Let |cc| be the norm of x e Z d , and |a"|* = 2 + |x|. 

Proposition 3.1 (Existence of stationary corrector and moment bounds [12]). 
Recall that we assume d ^ 3. For every A > 0, there exists a unique stationary solution 
to equation (3.2). Moreover, for every p^ 1, E{|</>a,{| p } and E{\D<f>\^\ p } are 
uniformly bounded in A > 0. The limit (f>£ = limA->-o is well-defined in L P (Q) and 
is the unique centered stationary solution to (3.3). 

Denote by G\(x,y) the Green function of A + V*a(x)V (the dependence on the 
randomness ( is kept implicit) and recall that G(x,y ) = Go(x,y). The following 
pointwise bound holds: 

for some c, C > 0. The following result controls the derivatives in the annealed sense. 
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Proposition 3.2 (annealed estimates on the gradients of the Green function [18]). 
For every 1 ^ p < oo, there exists C p < oo such that for every A ^ 0 and every e, e' e B, 

IvGa 

||VVG A (e,e')M^^. 

Remark 3.3. Notice that VG(x,e) (for x e Z d and e e B) denotes the gradient 
of G(x,-) evaluated at the edge e. Similarly, VVG(e,e') denotes the gradient of 
VG(-,e') evaluated at the edge e. 

3.4. Notation. We summarize and introduce some more notations used throughout 
the paper. 

• For i - 1,... ,d, e e M and /, Vi/(e) := Vi/(e). Recall that without any 
subscript, V/(e) = /(e) - /(e), and for x e Z d , Vif(x) = f(x + e*) - f(x) 
and V*/(t) = f (x-e z ) - f(x). 

• We write a < b when a ^ Cb for some constant G independent of e, e, a;. 

• For a, 6, c > 0, we write a < if for any 6 > 0, there exists Cg > 0 such that 
a ^ Cg pr j . In this way we have 

logM* < _J_ 

kls ~ Mr ’ 

• The Laplacian on Z d and the horizontal Laplacian on the probability space 
are both denoted by A = -V*V and A = -D*D. 

• For a random environment / and edge e € B, we obtain the environment 
perturbed at e by replacing ( e with an independent copy without changing 
other components (C e ')e'*e- The resulting new environment is denoted by 

C e - 

• For a random variable / and an edge e 6 B, the variable perturbed at e is 
denoted by f e (() := /(C e )- F° r a stationary random field f(x) = f(r x (), 
the field perturbed at e is denoted by f e (x) := f(r x ( e ). 

• The discrete homogenized Green function of A + V*ahV is denoted by 
Gh,A (x,y) for A ^ 0, and G^(x,y) = Gh.o(x,y). Recall that the hetero¬ 
geneous Green function of A + V*aV is denoted by G\(x,y ), and that 
G(x,y) - Go(x,y). The continuous homogenized Green function of -V-ahV 
is Gu(x,y). 

• {ei,i = l,...,d) represents the canonical basis of Z d , the corresponding 
edges, and the column vectors so that the identity matrix Id = [ei,..., e^]. 

• For functions of two variables, e.g., G(x,y ) with x,y e Z d , we use V x ,i, V y y 
to denote the derivative with respect to Xi , yi respectively. 

• The arrow => stands for convergence in law, and N(0,a 2 ) is the Gaussian 
law with mean 0 and variance a 2 . 

• a v b = max(a, b) and a a b = min(a, b). 

4. HELFFER-SJOSTRAND REPRESENTATION, TWO-SCALE EXPANSION OF THE 

Green function, and second order Poincare inequality 
We divide the proof of Proposition 2.5 into two steps. First, we show that 
(4.1) e- d Var {U< f) (g)}-*<T 2 g , 
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with a 2 defined in (2.3). If a 2 = 0, we conclude fyp\g) -*■ 0 in L 2 (Q). Next we 
assume cr 2 > 0 and show 

,4.2) „ 

VvmfwFo)) 

Once this is done, we can write 

% u \g) = u c f) (9)-^{Ue f) (g )} x yj e -d Ya , v { U (f\ g )} 

to conclude that (g) =>■ N(0,a 2 ). 

The proof of (4.1) uses the Helffer-Sjostrand representation and a two-scale 
expansion of the Green function, while the proof of (4.2) relies the second order 
Poincare inequality developed by Chatterjee [5] aud revisited in [22]. Both of them 
require taking vertical derivatives of lie ; ( g ) with respect to the underlying Gaussian 
variables £ e . Recall that U G \g) = e d+2 T. x , y ^z d G(x,y)g(ex)f(ey) is a finite linear 
combination of G(x,y). By Lemma A.l, we have 

d e U G) (g) =-s d+2 E d e a e 'VG(x, e)'VG(y, e)g{ex)f{ey). 

x,yeZ d 

We introduce the key elements in proving (4.1) and (4.2) in the following section. 

4.1. Helffer-Sjostrand representation and a two-scale expansion of the 
Green function. 

Proposition 4.1 (Helffer-Sjostrand representation [23]). Let /, <7 : -*■ R 6 e centered 

square-integrable functions such that for every eel, d e f, d e g € L 2 (fl). We have 

(f,9) = E(<9 e /,( 1 + -§0~ 1 d e 5). 

eeB 

Moreover, for every p ^ 2, (1 + is a contraction from L P (Q) to L p ($l). 

Since 

Ue f) (g) = £d+2 E G ( x ,y)9( £x )f( £ y), 

x,yeZ d 

the proof of (4.1) is reduced to asymptotics of Cov{G(x,y),G(z,w)} when the 
mutual distances between x,y,z and w are large. By applying Proposition 4.1 and 
Lemma A.l, the covariance is given by 

Cav{G(x,y),G(z,w)} = Y,(deG(x,y),(l + J^) _ 1 S e G(;r, w)) 

eeB 

= Y,{d e a e \7G(x, e)vG(t/, e), (1 + i c )^ 1 9 e a e VG(x, e) VG(w>,e)). 

eeB 

To prove the asymptotics, we need an expansion of X/G(x,e). The following proposi¬ 
tion is our main result in this section and one of the main ingredients to prove (4.1). 

Proposition 4.2. Recall that G and G h are the Green functions of V*aV and 
V*OhV respectively. For any eeB, we have 

|| VG(0,e) - VG h (e) - £ V fc G h (e)V<Me )|| 2 < l 0 p]§^. 
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An immediate consequence is that for any X € L 2 (H), 

(4.3) \(X, VG( 0, e)) - f V fe G h (e) (X, (e k + V&)(e))| < jXj 2 . 

fc =i la* 

By translation invariance of the environment, we further obtain for any x 6 

(4.4) [(A, vG(i, e)) - f V fc G h (e - s)(Jr, (e fc + V<k)(e))l < |X|| 2 b f lg ~ . 

fc=i Itl *|» 

Remark 4.3. (4.3) is an improvement of [23, Theorem 5.1]. 

To prove Proposition 4.2, we introduce the flux corrector, following [10]. For every 
i = 1,..., d, define qt = a{ei + D(fi) - a^ei, which describes the current correction, 
and qtj to be its j —th component. By the corrector equation (3.3), we have D*qi = 0, 
and by the expression of a\, in (3.4), E{ q Kj } = 0. 

We need the following integrability property. 

Lemma 4.4. Fix any i,j,k = 1,..., d. For A > 0, let o\ solve 

(A - A)ct a = D k q ij} 

then a\ is bounded in L 4 (i 2) uniformly in X,i,j,k. Furthermore, <j\ converges in 
L 2 (H) with the limit a e L 4 (Q) and solving 

-A a = D k qij. 


Proof. We first apply the spectral gap inequality in the form given by Lemma A.2 
to <j\ and obtain 

VeeB / 

Then we compute a\ - a e x for fixed e e B. Let Ga,a be the Green function of A - A, 
we have 

oa = cG(0)= Y, G A ,\(0,y)Vkqij(y), 

yeZ d 

°a = 5 a( 0 )= Z Ga,a(0 ,y)X7 k qtj(y)- 

yeZ d 

Since c/, 7 is the j- th component of a(ei + Difii) - ethe^ and cih = aid, we have qij = 
ajli = j + ajDjtpi - ali = j, which implies 

I Qijiy) - Qtj(y)\ $ l v=e{± +1 VjMy)\) + \VjMy) - I- 

Now we have 

kA - < Z |VfeG A ,A(0,2/)| (ly=e(l + \V j4>i(y)\) + I VjMy) ~ l) 

yeZ d 

=|VjG A ,A(0,e)|(l + |V^i(e)|) + Z |V£G a ,a(0, vWMv) ~ V,-# (y)\ 


so ^/E{|<ta - cr^| 4 } < ^/E{|/i| 4 } + \/E{|/ 2 | 4 }. The homogeneous Green function 
satisfies |V£Ga,a( 0,a:)| < \x\\~ d , thus \/E{|/i| 4 } < |e| 2_2d by Proposition 3.1. For / 2 , 
we write 

\h\ A = Z niVfeGA,A(0,2/„)||V,ki(2/n)-V j ^ e (y„)|, 

Vl,V2,V3,Vi^ d n=1 
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and by Lemmas A.3 and A.6 we have 


E{|/ 2 | 4 } < 



\ 4 

—- - — \ < - - -, 

hit 1 \y~e\*J |e|i 4d_4)_ ' 


so \/E{|/ 2 | 4 } < l/|e|i 2d . In summary, we have 

1 


E K) s E . |(M - a) . 

VeeB |e|* 


and since dZ 3, we conclude E{cr 4 } < 1. 

To show the convergence of a\ in L 2 (fl), we only need to prove that (cta 1 ,(Ja 2 } 
converges as Ai,A 2 -»■ 0. By the Green function representation, we have 


(<L\ !,o-a 2 >= E ^kG AM (0,y 1 )VtG AM (0,y 2 )E{q lJ (y 1 )q lJ (y 2 )}. 

yi,yaez d 


By Lemma A.4, \E{q i j(y 1 )q ij (y 2 )}\ < ]yi _l 2] d- • Furthermore |V£G AjA ( 0, y)| < |y|i d , 
by the dominated convergence theorem, we have 

(oAi,0A 2 )^ E VfeG A ,o(0,i/i)VfcGA,o(0,t/2)E{g ij (yi)5y(i/2)}- 

2/i,y 2 eZ d 

Therefore, <j\ converges in L 2 (£i). Its limit a is in L 4 (fl) by Fatou’s lemma. By 
sending A -*■ 0 in (A - A)a\ = D^qij, we obtain -Act = Dj-q^, and the proof is 
complete. □ 


We can now define the flux corrector {<Jij k ,i,j,k = 1,... ,d}: 

Lemma 4.5. There exists a tensor field {&ij k ,i,j,k = l,...,d) such that 

* &ijk ~ 

* &ijk ^ (^); 

* = D kqi j — D jQik and ^£>=1 D^<7ij k = Qij. 


Proof. For every i,j,k = 1,..., d and A > 0, we consider the equation 

( 4 -5) (A - A )af jk = D k qij - Djq lk . 

Lemma 4.4 ensures that E{|erA fc | 4 } < 1, that crA fc converges in L 2 (fl), and denoting 
the limit by <Jijk, we have <Jij k e L 4 (fl) with -A<Jij k = D k qij - Djqn.. 

The skew symmetry o t j k = -<r lkl is clear by (4.5). 

To show Efc=i D k (iijk = qij, it suffices to prove A(Efc=i Dl<Jij k - Qij) = 0. Indeed, 
D(T,k=l D k Cijk ~ Qij ) = 0 implies D k &ijk ~ Qij ~ const by ergodicity, and since 
E {£fc=i D* k a ijk ] = E {q i:j } = 0, we have '£, k =i D* k a ijk = q t j. Now we consider 

A ( E D k a ijk - Qij ) = lim A ( E D *k°ijk ~ Qij) 

k= 1 A_>U fc=l 

d 

= I™ E D k( D jQik ~ D k qij + \a* jk ) - Aq lr 

A ^ U fc=l 

Since Y,t=i D k qi k = 0 and cjA fc is uniformly bounded in L 4 (fi), we have A(2 )Ei -D k a ijk~ 
q^) = 0, and this completes the proof. □ 

Proof of Proposition 4-2. We follow the proof of [23, Theorem 5.1], but use the flux 
corrector as in [10] to simplify calculations. Define 

z(x) := G(0,x) - G h (x) - Yj V/cGh(x)^ fc (x), 

k= 1 
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as the remainder in the two-scale expansion of the Green function, the matrix 
function R:= -[qi,... ,q d ] by 

f / ( d ') — ~Qji( x ') — ^1 i=j ~ C/('T ) ( ^'{ y A Vj0j(x)), 

and h : Z d -»■ M by 

d ( d 

h(x) = - XI a*( x ) E &i( x + Zi)ViVjGu(x) 

i= 1 \ i=l 

By [23, Proposition 5.6], we have 
(4.6) z(x) = E G ( x ,y) E Rijiy ~ ei)SJ*SJjG h (y) + E G(x,y)h(y). 

ye Z d i,j =1 yeZ d 

Consider the first term on the right-hand side of (4.6). Since R l3 - -q 3 i - 
- Efc=i Dk a jik by Lemma 4.5, we can write 
d d 

E R ij (y-e i )V*V j G h (y)=- E Vfe(y - e i )V 1 *V J -G h (y) 

i,.7=l i, < 7,fe=l 

d 

= - E Vfc(CTj ifc (y-ei)Vi VjGh(y)), 

where the last equality uses the fact V* k (f{x)g( cc)) = \I* k f{x)g{x) + f{x-ek)V* k g{x) 
and + erjfci = 0. Therefore, by using the flux corrector a, we can write 
T,i,j=i Rijiu ~ e i) Vj" \7jGb(y) in divergence form. Note that h is in divergence 
form. An integration by parts leads to 

d d 

z(x)=- E E Vy,fcG(x,y) E Vjikiy ~ c)V* VjGh(y) 

yeZ d &=1 *,j=l 

d d ~ 

- E E v y4 G ( a; >2/)^(y) E^(y + e i)v i v J G h (?/), 

2/6 Z d 1=1 J = 1 

so for e e B, we have 

d d 

Vz(e)=- E Y VS7 y,k G ( e >y) E d jik(y - e i)Vi VjG h (y) 

yeZ d fc=l i,j=l 

d d __ 

_ E E vv y.i G ( e ’2/)«i(y) E^(^ + e i)V l V J G h (y). 

i/6Z d 1=1 j=i 

Note that 


V; 2 :(:r) =V;G(0,x) - V;G h (:r) - ^ (V^/ c (x)VfcG h (x) + + e*)ViV/cG h (x)) 


fc=i 


= ( ViG(0,x)-ViGh(x)- ^ V^fc(x)V fc Gh(x) I - ^ 4> k (x + ei)ViVfcG h (x), 


fc=i 


)-E 


and moreover, by the moments bounds on tf>]. provided by Proposition 3.1 and the 
fact that |ViVjGh(ir)| < |a;|* , we have 

||^(* + e i )V i V fc G h (a :)||2 < r4 < 

\x\i \x\i 

As a consequence, in order to prove Proposition 4.2, it is enough to show that 


(4.7) 


|Vz(e)|| 2 < 


log |e[. 
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Ill order to prove (4.7), we note that Vz(e) is a finite linear combination of terms in 
the form £ y6Z d VV y , fe G(e, y)f(y) V ■ VjG h (y) or VV y , fc G(e, y)f(y) V; VjG h (y) 
for some i,j, k and /. Clearly, they can be bounded by D ys z d |VV y ,fcG(e, y) f (y)\\y\* d , 
so we have 

|v*(e)|| 2 <| £ |VV y ,feG(e, y)f(y)\\y\* d \\2 

yeZ d 

< E lyl* d |VV y ,/cG(e, y)f(y)\\2 

ye Z d 

< E \y\* d \\Wy,kG(e,y)\\ 4 \\f(y)U. 

ye Z d 


When / = djik or di(j)j, ||/(y )||4 is uniformly bounded by Lemma 4.5 and Proposi¬ 
tion 3.1, thus by applying Proposition 3.2 and Lemma A.6, we obtain 


|V*(e)|| 2 < £ 

yeZ d 


1 1 

\y\* \y-§$ 


< 


log [eU 

|e|* 


The proof of Proposition 4.2 is complete. 


□ 


4.2. Second-order Poincare inequality. Let cLk be the Kantorovich-Wasserstein 
distance 

d K (X,Y) = sup{E{/i(X)} - E{h(Y)} : sc 1). 

In order to show that the rescaled fluctuations are asymptotically Gaussian, we will 
use the following result. 

Proposition 4.6. [22, Proposition 2.1] Let F e L 2 (Ll) be such that E{F} = 0 and 
E{F 2 } = 1. Assume also that F has weak derivatives satisfying Y, e E{|9 e F| 4 } 5 < 00 
and E{\d e d e 'F\ 4 } < oo for all e,e' eM. LetY~N( 0,1). Then 


(4.8) 


d K (F,Y)<\/~^ 



||d e F||4||d e d e ,F||4 - 


Using the above result, we only need to show the following lemma to prove (4.2). 

Lemma 4.7. Let 


(4.9) 


« 2 := 


E (iWPitiU E \9ed*uP(g)u) ■ 

e'eB V eeB / 


If dg defined in (2.3) is not zero, then 


Yar 2 {U^\g)} 


< e d \ log er| 2 . 


Remark 4.8. By Proposition 4.6 and Lemma 4.7, if a 2 t 0, then we actually obtain 
the convergence rate 


d K 


U { J\g)-E{U { J\g)} 
^Var {U { J\g)} 


\ 


,N( 0,1) 


< e 2 1 log e|. 


/ 
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5. Convergence of the variance 


The aim of this section is to prove (4.1). 

Recall that ui J \g) = e d+2 Y. x , y ez d G ( x ,y)g( £x )f( £ y), so 

Yar {U G) (g)} =e 2d+A Z Cov{G(x, y), G(z, w)}g{ex)f{ey)g{ez)f{ew) 


=£' 


x,y ,z,w€£j 

2d+4 V- y \ W > 


Z Cov{G( —, -), G(-, — )}g(x)f(y)g(z)f(w). 


x,y,z,w€e TL d 


£ £ £ £ 


The covariance is given explicitly by the Helffer-Sjostrand representation 
(5.1) Cov{G(a:,y),G( 2 ,u;)} = Y,( d e G (x,y), (1 + _S?) _1 d e G( 2 !, w)), 


ee® 


and since d e G{x,y ) = -d e a e VG(x, e) VG(y, e) by Lemma A.l, (5.1) is rewritten as 
Cov{G(a:, 2 /),G( 2 ,iy)} 

(5.2) = ^( 9 eCle vG(x, e) VG(y, e), (1 + &)~ 1 d e a e yG(z, e)VG(w, e)). 


ee® 


To prove the convergence of £ d Var{£t£ (<?)}, we use the two-scale expansion of the 
Green function obtained in Proposition 4.2. For e e B, Xi e Z d , i = 1,2,3,4, define 

1 


(5.3) 

and 




ee®i=l I— ^i\* j=l,j±i I— *^jl* 


e- x- 


Id-1 


K.(x,y,z,w) := Z Kijki Z ViGh(u-x)VjGh(i'-y)VfeGh(w- 2 :)VzGh('!;-ui) 

i,j,k,l =1 veZ d 

with K ijk i given by (3.6). 

Proposition 5.1. |Cov{G(ic,y),G(z, w)} -lC(x,y,z ,«;)| < £(x,y,z,w). 

Proof. Each term in (5.2) contains four factors of gradient of the Green function. 
We first consider VG(i,e) and let 

A' = <9 e a e VG(y, e)( 1 + if) _1 9 f ,a e VG(z,e)VG(w,e). 

By (4.4), we have 

|(A, VG(z, e))| - Z VfcGh(e - x)(X, (e fc + V«(e))| < ||X|| 2 l0g|g ~ ^ . 

fc=i |e-x|; 

By Proposition 3.2 and the fact that (1 + Jz?) -1 is a contraction from L P (Q) to 
L p (£l) for any p > 2, we have \X\ 2 <\e-y\\~ d \e-z\\~ d \e-w\\~ d , so 

d 

|Co v{G(x,y),G(z,w)} - Z Z VfcG h (e - x)(X, (e fc + V</>fc)( e ))l < £(x,y,z,w). 

ee® fc=l 

Now we carry out the same argument for VG(j/,e),VG(z,e),VG(w,e), and in the 
end obtain 


|Cov{G(a;, y),G{z, w)} - IC(x,y,z,w)\ < £(x,y,z,w). 
The proof is complete. 

Proposition 5.1 leads to 

|£- d Var{z4 /} ( 5 )} ~ d+4 


□ 


£“* Z /C(-, -)g( x )f(y)g(z)f(w)\ 
£ £ £ £ 


x,y ,z,w££ lA 

<e d+4 Z £(-,-,-,-MaO/(yM^)/(w)|. 
£ £ £ £ 

X,y,Z,W££ 
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Hence, the proof of (4.1) will be complete once we have proved the following two 
lemmas. 

Lemma 5.2. e d+i Y, x , y>ZtWee zd £(*, f, §, j)\g(x)f(y)g(z)f(w)\ -> 0 as e - 0. 

Lemma 5.3. e d+i 'L x , y , z , w ^z d /C(f, f, §, ™)g(x)f(y)g(z)f(w) a 2 as e ^ 0. 

In the following, we assume |<?|, ]/| < h for some h e C c (R d ). 

Proof of Lemma 5.2. By denoting the number of different elements in {x,y,z,w} 
by s, we decompose Y,x,y,z,weeZ* = £ l«=i + E l s =2 + E R=3 + E 1«=4- The following 
estimates are obtained with an application of Lemma A. 10. 

When s = 1, 

e d+4 IW:(-A £ ,-)| g(x).f(y)g(z)f(w)\<e d+4 X h(x) 4 . 

£ £ £ £ 

Since h e C c (R d ), e d 'Z xie z d ^(*) 4 -*■ / R d h(x) 4 dx which is bounded, so we have the 
r.h.s. of the above display goes to zero as e -> 0. 

When s = 2, 

,x y z w , 


e d+4 X 1 s= 2 £(~, -)lff(^)/(?/)5(^)/(M’)l 

£ £ £ £ 

\2\ X V\2-2d 


x y 11—rf 


~ d+4 X (^) 2 My) 2 + ^) 3 My))l---li' d 


<, d+4 X K x YKy) 2 \-- -\l~ 2d + M*) 3 Mj/)I- - - 

x*y6eZ d £ £ £ £ 

For x *0 £ Z d , we have |x|* > |x| ^ 1, so the r.h.s. of the above display is bounded by 
£ d 

x±yzeZ d 

^£ 2d+3 Y (h(x) 2 h(y) 2 + h(x) 3 h(y))\x-y\ 1 - d . 

xtyeeZ d 

Similarly, e 2d ff x ± y <ieZ d {h{ x ) 2 h{y) 2 + h{x) 3 h{y))\x - y| 1_d converges as a Riemann 
sum, which implies 


_2d+3 


Y (HxfKy) 2+ K x YKy))\x-y\ 


1 -d 


0 


x±yeeZ d 

as e -»■ 0. 

The discussion for s = 3,4 is similar to s = 2, so we omit the details. 
When s = 3, 

_d+4V, cr x y z w ' 


e XI 1 s=3£(~, -)IS'(a:)/(y)fl , ('2)/(^)l 

£ £ £ £ 


<£' 


3d+2 


Y h( x ) 2 h(y)h(z) I 


l 


x,y,zes Z“ 


y| 


d-l 


1 


1 


X - X 


d-l 


Iz/- 


• I d-l 


' £ 2 . 


|{®i2/i«}I=3 
When s = 4, 

e d+4 X !s= 4£(-, -> -)l5(^)/(y)5(^)/(^)l 

£ £ £ £ 

< e (4d+i)- £ l s=A h(x)h(y)h(z)h(w)~ 


i- 


-y|(d— 1)- 


- 1 )— 


;|(d-l)- 


The proof is complete. 

Proof of Lemma 5. 3. Recall that 


□ 


fc(x,y,z,w) = X K ijki X ViG h (u-x)VjGh(?;-y)VfcGh(u-2)VzGh(ii-iu). 
i,j,k,l=l ueZ d 
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By defining T ijk i(v,x,y, z,w) := VjG h (t; - x)S7jG h (v - y)VfeG h (« - z) V;G h (v - w), 
we only need to show the convergence of 

_d+ 4 £ XI Fijkiiv,)g(x)f(y)g(z)f(w) 


*ijkl 


= £ 


£ £ £ £ 


x,y,z,w££Z d veZ d 

for fixed i,j. k, l. 

We claim that VjG h (f- f ) can be replaced by d Xi Gb(v- f) in Tijki{v, f, f 7 ) 
of the above expression with the sum over w =£ x/e. Indeed, by [23, Proposition A.3], 
for x + 0 , 

|VjGh(x) - d Xi Q h (x)\ < \x\~ d . 

If we define J^ kl (v,x,y,z,w) := d x , Gb(v - x)V,G h (r - 2 /)VfcG h (u - z)V/G h (u - u>), 
the error induced by the replacement can be estimated as 

IA?M- £<i+4 E E Fju( v , -> -)g(x)f(y)g(z)f(w)\ s -A + J 2 , 

x,y,z,w£EZ d v±j £ £ £ £ 


with 


and 


Ji = e' 


d+4 


y > . / X X Z 'UJ \ \ , / \ 1 / \ 1 / \ 1 / \ 

E I Fijkl{~, — )l h{x)h{y)h{z)h(w). 

x,y,z,weel, d £ £ £ £ £ 


x,_ d \ v _y\i- dl 
£ £ 


il-di 


J 2 = e d+4 e e i«- -r> - -ii - > - -irv- -1 i- d h(x)h(y)h(z)h(w). 

— *-“* cr C C C 

x,y,z,weeZ d vt- 

For Ji, by using |VjGh(x)| < |x|*~ d and considering different cases according to 
whether y,z,w = x as in the proof of Lemma 5.2, we obtain 


1 -d. 


Ji<e d+4 Z 

x,y,z,we 

For J 2 , we note that T, v *^ I 11- ^\ d \ v ~ f 
£(x,y,z,w) defined in (5.3), so we can apply Lemma 5.2 to show J 2 -*■ 0. The claim 
is proved. 

By following the same argument for V 7 Gh(w - -), VfcGh(u - -), V;Gh(u - —), we 
derive 

I hjki-e d+4 E E r jkl (v,-,-,-,-)g(x)f(y)g(z)f(w) |-o, 

x,y,z,wzeZ. d v*-, 11 -,-,— £ £ £ £ 

’ ’ e ’ e ’ e ’ e 

with T ljkl (v,x,y,z,w) := d Xi Gb(v - x)d Xj Gb(v - y)d Xk Gb(v - z)d Xl Gb(v - w). Since 
dZ 3, Gb(x) = Ch]x| 2_d for some constant Ch, and d Xi G h(x) = Ch(2 - d)x i /|x| d , so we 
have 


,x- y a - d ,x- z a _ d .x - w a _ d 


£ £ 

x \-d\ n . _ y\l-d 


h(x)h(y)h(z)h(w) -*■ 0 . 


v-^\ 1 - d \v--\\~ d ^ with 


^d+4 


£"■’ E E r^ k \v,- y- Z -^)g{x)f{y)g{z)f{ W ) 

x,y,z,w££ £ £ £ £ 

=£ 5d Z 1 v*x,y,z,wF‘ jkl (v,x,y,z,w)g(x)f(y)g(z)f(w) 
x,y,z,w,ve£ Z d 

/ T ljkl {v,x 1 y,z,w)g(x)f(y)g(z)f(w)dxdydzdwdv. 

Js? d 

The proof is complete. 

6. Convergence to a Gaussian when er 2 > 0 
Recall that in order to prove (4.2), that is, 

U { J\g)-HW\g)} ^ 

\J\ax{Uy\g)} 

we only need to show Lemma 4.7. 


□ 


N{ 0,1), 
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Proof of Lemma f. 7. We first prepare the ground by estimating the terms appearing 
in the definition of re 2 . By a direct calculation, we have 

d e U { J\g) = -e d+2 Y d e a e VG(x,e)VG(y,e)f(sy)g(£x), 

x,yeZ d 

and 

d e 'd e U( J) (g) =-e d+2 Y dla e VG(x,e)\?G(y,e)f(£y)g(£x)l e ' =e 

x,yeZ d 

+ £ d+2 d e a e d e 'a e 'VG(x, e 7 ) VVG(e, e')VG(y, e)f(£y)g(£x) 

x,yeZ d 

+ £ d+2 Y d e a e d e 'a e 'VG(x,e)VG(y,e')VVG(e,e')f(£y)g(£x). 

x,yeZ d 


By Proposition 3.2 and the fact that a e = rj(Ce) with |t? / |, |'?7 /, | uniformly bounded, 
we have for any p ^ 1 (with the multiplicative constant depending on p ): 

\\deW\g)\\p <£ d+2 Y ] -- l ~uzi\f(zy)g(zx)l 

x.vOA F-e|S \y-e\i 


and 


\\de'd e uP(g)\\ p <£ d+2 Y 7---wrrl/( £ 2/)5(ea;)|l e '=e 

x.v& f-lI* \y-Q* 

£ \ x - |e - §L\ * I y - g\t x 

L\ x -^\e-^\ y-lf\^ f{£y)9i£X)l 


x,yeZ d 
d +2 


x,yeZ d 

Since f,g are both bounded and compactly supported, we apply Lemma A. 11 to 
obtain 

< i 12 

( 6 . 1 ) 

and 


\\d e U^\g)\p < £ (£ , _ | d _i 1 

\xti. d \ x LI* / 


ee 


\2d-2 ’ 


( 6 . 2 ) 


\d e ’d e U^\g)\ p <£ (£ | _ | d _r l|®|se-i] ( le '= e + ip _ e /id) 

Wz d F LI* / IL <LI* 


We are now ready to estimate re 2 . By (6.1) and (6.2), we have 
« 2 = £ (£ \\d e ’d e U { e f) (g)U\\d e U { e s) (g)\\ 

e'eB VeeB / 

l 



< 

Applying Lemma A. 11, we get 

_2 d 


\e-d\ d |ee| 2d_2 



——) < y 

-eL\ d J ha 


£ id \ log e|' 


\ee 


r\2d 


^ 4d |iog £ | 2 £ 


Xil d + l e:r l) 


2d 


<e 3d | loge|' 
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To sum up, k 2 < £ 3d |loge| 2 . By (4.1), Var {Ue f \g)} > e d if a 2 t 0, which leads to 

k 2 e 3d ]loge| 2 dl ,2 

< — ' _ 6 ' < e d \ loge| 2 -+ 0, 


Var 2 {G/H (g)} 
and the proof is complete. 


r 2 d 


□ 


7. Tightness in C 1 “ c 

Roughly speaking, for a < 0, a distribution F is a-Holder regular around the 
point x 6 M. d if for every smooth, compactly supported test function %, we have 

(7.1) F[e- d X (e~\--x))]<e~ a (s-0). 

By [21, Theorem 2.25], in order to prove that ( $/^' > is tight in Cj“ c for every a < 1- 
we only need to prove the following proposition. 

Proposition 7.1. For any g e C c (R d ), let g\(x) = \~ d g(x/ A) . For all p ^ 1, there 
exists a constant C = C(p,g ) such that 

(dx)\\p $ 

Proof. We follow the proof of [22, Proposition 3.1] . Recall that (g) = \g) - 

E{U^ f \g)}) : by writing 

^' £ (/) (5a)a^ _ 1 = L-j- Z Z ( G ( x >y)- E { G ( x >y)})f( £ y)9(^)--x^x, 

A 2 x^L d yeZj d 

and we need to show that X e \ is uniformly bounded in L p (fl). Since E{X e >} = 0, 
in particular 

. 2 


(7.2) 




< 1 


holds for p = 1. We argue inductively, assuming that (7.2) holds for some p = n, and 
showing that it also holds for p = 2 n, which would complete the proof. 

Since E{X 2 "} = E{X £ " A } 2 + Var{AT” A } < 1 + Var{AT” A }, it suffices to show 
Var{X™ A } < E{X 2 "} 1- ". By the spectral gap inequality (see [23, Corollary 3.3]), 
we have 

Var{X" A ) 4 Z E{|5 e A- A | 2 } = Z E{| nX^d e X e , x \ 2 } 


eel 


eeB 


<E{|X £iA | 2li } 1 ^Z E {l^^,A| 2 "}"- 

eeB 

So we are left to prove £ e(E B E{|9 e A' eiA | 2n } « < 1. Since 


-+2 

f 2 fnr 

d e X etX = --j— Z 3 e fleVG(a;,e)VG(i/,e)/(£i/)s(y) 
A2 +i A 

i )d\ 


and /, g e C c (R d ), by applying Proposition 3.2 and Lemma A.11 we obtain 

dj 


\\d e X £ , 


A 2 n 


F 2 +2 

<— Z 

4 + 1 , 




X€l 


\ x - e \d-l ■*‘b|S7 Z 


1 


\y-e\ 


d -i ti/l*; 


£2 +i A 1 1 1 


'a! + 1 £ IfelZ 1 e |ee|f 


HZ 


If e|* _1 M* -1 ’ 


which implies 

EE{|a,x eiJ j 2 "}i < y (£)‘'_i 

eeB A ' I \ Ll* 


2d-2 \pp\2d-2 

A—I* l CC l* 


e(()‘ 

eeB ' A ' 


\Lp\2d-2 

I \Sl\* 


< I- 


The proof is complete. 


□ 
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Appendix A. Technical lemmas 

Lemma A.l (vertical derivative of G(x,y ) with respect to £ e ). For e € M,x,y 6 
Z d , uie.fl, we have 

(A.l) d e G(x,y ) = -d e a e VG(x,e)VG(y,e). 


Proof. Fix e,x,y,f. By definition, the Green function G(x,y) = / 0 °° qt{x,y) dt with 
the heat kernel q t (z,y ) solving the following parabolic problem 

d t qt(z,y ) = -V*d(z)Vq t (z,y),z e Z d , 

with initial condition qo(z, y) = 1 z=y . We take d e on both sides of the above equation 
to obtain 

d t d e q t (z,y ) = -V*a(z)Vd e q t (z, y) + d e a e Wq t (e, y)(l z= e ~ l z =e), 
with initial condition d e qo(z,y) = 0. So d e qt(x,y ) is given by 

d e q t (x,y)= f Y q t - s (x,z)d e a e Vq s (e,y)(l z=e -l z=g )ds 

- - d e a e / SJq t _ s {x,e)SJq s {e,y)ds. 

Jo 

This leads to 

J f' oo r~ OG r~t 

d e qt(x, y)dt = - d e a e / / \7q t - s (x,e)S7q s (e,y)dsdt 

o Jo Jo 

= - d e a e VG(i,e)V%e), 

where we used the symmetry q t (x,y) = qt(y,x) in the last step. The proof is 
complete. □ 

Lemma A.2 (Spectral Gap Inequality to control fourth moment). For any f with 

E{/} = 0, 

(A.2) E{/ 4 }<(^\/E{|/-m) ■ 

VeeB / 


Proof. By [8, Lemma 2], if we define E e {/} := E{/|{£ e /} e '* e }, then 

E{/ 4 }<E{(y|/-E e {/}] 2 ) 2 }. 

eeB 


By expanding the right-hand side, we obtain 

E{/ 4 }< I E{|/-E e {/}| 2 |/-E e ,{/}] 2 } 

e,e'eB 


$ I \/E{|/ - E e {/}| 4 }E{|/ - E e /{/}] 4 } 

e,e'eB 

= (E\/E{l/-E e {/}| 4 }) • 

VeeB / 


It thus suffices to show that E{|/ - E e {/}| 4 } ^ |E{|/ - / e | 4 }. In order to do so, we 
write 


/-r = /-E e {/}+E e {/}-r, 


and observe that 

E{|/ - ri 4 } =E{|/ - E e {/}| 4 } + E{|E e {/} - / e | 4 } + 6E{|/ - E e {/}| 2 |E e {/} - /f } 
+ 4E{(/ - E e {/})(E e {/} - n 3 } + 4E{(/ - E e {/}) 3 (E e {/} - / e )}. 
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By first averaging over ( e (resp. £'), we see that the third (resp. fourth) term on 
the right-hand side is equal to zero, so the proof is complete. □ 

Lemma A. 3 (Sensitivity of gradient of correctors with respect to ( e ). For e e B, x € 
Z d , i, j = 1,..., d and p ^ 1, we have 

(A.3) E{|Vj0i(a;) - | p } < \x - e\~ pd . 


Proof. By the convergence of Vj4> x ,i(x) -*■ Vj0i( x) in L p ( fl), we only need to show 

E{|V^ A ,i(a:) - VAU X )\ P } ~ \ x -e\* pd , 

where the implicit multiplicative constant is independent of A. 

We write the equation satisfied by and tfi xi as 

X4>x,i{x) + S/*a(x)\/(j)x,i{x) = -V*a(a:)e i , 

M\ t ii. x ) + S/*d e {x)S/^ e Xi {x) = -V*d e (x)ei. 

A straightforward calculation leads to 

<t>\,i( x ) ~ 4>%ii. x ) = E g a (x,y)V*(d e (y)-a(y))(Vf> e x Ay)+ei), 

ye Z d 

so we have 

Vj<A ;~ ^j4>\,ii. x ) = E V x>j S7 v G x (x,y)(d e (y) - d(y))(V4> e Xi {y) + ef). 

yeZ d 

Since a e (y ) - a(y) = 0 when y + e, we conclude 

IvAi(z) - V^A,i(a:)l S |V x ,iV a G A (a:,e)||V^A,i(e) + e*|. 

By Propositions 3.1 and 3.2, the proof is complete. □ 

Lemma A.4 (Covariance estimate of gy,-). For i,j = l,...,d and x e Z d , we have 
(A.4) | E{ ^(0)^ (a: )}|<^k. 


Remark A.5. Similar results in continuous setting are given in [16, Proposition 4.7], 
Proof. By [16, (4.4)], we have 

|E{%(0)g ii (a:)}| =|Cov{q li (0),% i (ai)}| 

eeB 

Recall that qij - ajli = j + ajD.jcfi - dl, =) , so for e e B,x € Z d , we have 

fe(z) - Qij(x) I S I dj(x) - a*(x )|(1 + \Vj(f>i(x)\) + \S7j4>i(x) - V^-(a;)|. 

By Proposition 3.1 and Lemma A.3, we have 


which implies 


^{\qij(x) ~ Qt,(x)\ 2 } Z U=e + \ 


\®{Qij(o)qij(x)}\ < E 


, 1-2 d 


log |a 


eel 


e “ x - e\\ 


where the last inequality comes from Lemma A.6. The proof is complete. 


□ 
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Lemma A.6 (Estimates on discrete convolutions). For a,/3> 0 with a + (3 > d, we 
have 

(A-5) Yj yTj 7g < C a ,pF a tf(x) 

y ^d |2/|* \x - y\f 

for some constant C a j( g > 0 and 

1 1 Jog |#| * 

(A.6) F a> p(x) = y d l a v/3<d + QA/3 lav/3><i + A/3 l a v0=d- 


Remark A. 7. We will usually replace by 

Proof. The proof is standard. Since a + f3 > d, we only need to consider the region 
|x| > 100. 

For fixed x, let h = {y ■ ]?/| \x\/2}, I 2 = {y ■ \y-x\ ^ \x\/2}, and I 3 = Z d \(I 1 uI 2 ). 

We control the sum in each region separately. The proof for each case is similar and 
we only use the following two facts: 

• |x|* < \y - x|* in I 1 , |x|* < |y|* in I 2 and \y - x\* ~ |y|* in I 3 , 

• for any 7 > 0, E| y |<| x | \v\ * 7 ^ \ x \* ” 7 I 7 <d + log |x|* l 7 =d + l 7 >d- 
In I 3 , we have 


Y— _L^< v 

i 3 \v\?\x-y\% | wWx | /2 


-I~ lx|“ +/3 “ d ' 


If a v f3 < d, the discussion for /1 and I 2 are the same. Take /1 for example, we 
have 

Y— _I_< _L V _L<_I_ 

it \y\* \ x - y\* ~ M* M^|/2 ll/l? ~ |a:|* +/3_<i 

If a v (3 > d, we assume a ^ (3, so a > d. By similar discussion, in I 2 we get the 
estimate 


v 1 1 „ 1 1 , !ogl*k . - 1 

\ Mi \ x - yf. ~ \x\r s - d f<d w; e ' d e,d ' 


In 7i, we have Z|y|^|x |/2 M* a |# _ 2/1* ~ |x|* . So an overall bound is given by 


\x\/ = ld* QA/3 . 


If av/3 = d, we assume a = d. If f3 = d, in both /i,/ 2 we get the bound |a:|“ d log|a;|*. 
If f3 < d, in Ii we get a bound \x\ff log|x|* and in I 2 we get the bound |x|* _a_/3 = \x\f^. 
Therefore, an overall bound is \x\~^ log|x|* = |x|*“ A/3 log |m|*. 

The proof is complete. □ 

Lemma A. 8. Let xi, ... , x^ e 7L d be mutually different, and for each i, let i be such 
that \xj - Xi\ = min,*, \xj - xf\. Assume a,\, ..., a& e (0, d) and a* + aj > d for i + j, 


e n^-E^En |T W,, ‘ , 

ieZ d i d* i—lj±i i I j d* \%i 


Proof. For each i, we define the region I, = {y ■ \y - Xi\ < min 7 | y - Xj |}, i.e., the set 
of points that are closest to Xi. If y e we have | y - x 3 1 ^ |xj - Xj\/2 for any j ± i. 
Therefore, 

J?, y \y - XjW ~ y, \xj - Xi]^ d |y - Xi\f d I y - xftf 1 ' 
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Since + oq > d , the sum over y can be bounded using Lemma A.6, e.g., when 
cti < d for all i, we have 


in 

v*h i 


I y- 


n 


u 3\* 


j±i,i \ X 3 ' 


i+otj—d ' 


The proof is complete. 


□ 


Remark A.9. From the proof of Lemma A.8, we see that the condition on + aj > d 
for all i + j is not necessary to obtain similar estimates. For example, for each i, as 
long as we can find j t i such that on + aj > d , the integral in R can be controlled 
by a similar bound. 


Recall that the error function £ in Proposition 5.1 is given by 

1 


c( \ v 1 V 1 1°§I- Xi \* 

£{X 1 ,X2,X 3 ,X i ) = 2^ L 


n 


(A.8) 


‘ * ‘ * I p _ on . I d ^ I s} _ rn . I d 1 

eel 2=1 \- j=l f j±i I- ' l 3\* 


II 


1 


l 2=1 


4 

w n 

■I* 3=1,3*^ 


1 


\V ~ X 


I <2-1 ‘ 
3 1* 


By using Lemmas A.6 and A.8, we have the following control on the error function: 


Lemma A. 10 (Estimation of £(x, y, z, w)). Let x,y, z,w € Z d , 

• if x = y = z = w, £(x,y, z,w) < 1. 

• if x = y = z * w, £(x, y, z, w) < \x - w\ff d . 

• if x- y + z- w, £(x, y, z, w) < \x - w\ff 2d . 

• if x = y and y,z,w are mutually different, let S = {x,z,w}, 


£(x,y,z,w) <£ n 

vzS ueS\{v} 


1 

\u-v\ff 1 


• if x,y,z,w are mutually different, let S = {x,y,z,w}, 


£(x,y,z,w) < I n 


,(d-l) ' 

veSueS\{v} \U- U|i 


Proof. The proofs of different cases are similar. We only discuss the case when 
x,y,z,w are mutually different. By (A.8), we consider a term of the form: 

y 1 1 1 _ 1 

v iz d \ v ~ x \*~ ~ y I* -1 1^ - ~l* _1 \ v - w \*~ 1 


Recall that S = {x, y, z, w}. For u 6 S, let I u = {v e Z d : \v - u\ = min^s \v - g|}. The 
proof is the same as in Lemma A.8, so we do not provide all details. 

For I x , we have 


I 

vel x 


\v - X “ \v ■ 


y IS 


tI 


\x - z\ d 1 |a; - W|“ * \v - a; is \v ■ 


y\ 


< 2-1 


n 


, |(< 2 - 1 )-' 
utS ^{ x } \X ~ U \\ 


Note that we replaced 


\x-z\ff 1 \x- 


d=r by 
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For I y , we have 


1 


1 


1 


1 


y 

v7i y \ v - Ai K’ - y\i 1 \v - Ai 1 k - Mi 1 

:_I_I_y_I_I_ 

' I y - At 1 1 y - ML 1 Jkd \v - x\i~ \v - ylf 1 


n 


i 


ueS*{y} \y-u\{ d 1} 

The sums in I Z ,I W are discussed in the same way. The proof is complete. 
Lemma A.11. For x € Z d , p > 0 and e e (0,1), 
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We refer to [22, Lemmas 2.2 and 2.3] for a proof. 
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